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Abstract. We define optimal per-particle fluctuation and correlation measures, relate 
fluctuations and correlations through an integral equation and show how to invert that equation 
to obtain precise autocorrelations from fluctuation scale dependence. We test the precision of 
the inversion with Monte Carlo data and compare autocorrelations to conditional distributions 
conventionally used to study high-pt jet structure. 



1. Introduction 

Fluctuations in nuclear collisions measured at a single bin size or scale could arise from many 
possible configurations of a multiparticle momentum distribution and are therefore difficult 
to interpret. However, the scale dependence of fluctuations over a significant scale interval 
does contain detailed information about multiparticle correlations which can be extracted with 
the proper techniques [Q. Information about the absolute location of event- wise structure 
is lost, but those aspects depending only on position difference are retained in the form of 
autocorrelation distributions 0. In this paper we consider fluctuations on binned momentum 
space, Pearson's correlation coefficient and autocorrelation density ratios for multiplicity n and 
transverse momentum pt. Autocorrelations can be inferred directly from pair ratios or from a 
fluctuation/autocorrelation integral equation which we derive. Inverting the integral equation we 
obtain autocorrelations which can be compared with more conventional conditional distributions. 



2. Fluctuations and correlations on binned spaces 

Correlation analysis of nuclear collisions reveals information arising from event-wise changes in 
multiparticle momentum distributions. The data system is an ensemble of event-wise particle 
distributions on momentum space (pt, rj, 4>) or (yt, rj, (j>), where pt is transverse momentum, wit is 
transverse mass, w is pseudorapidity, <f) is azimuth and yt = ln{(m^ + Pt)/'mo} is transverse 
rapidity with pion mass assigned to niQ. The momentum space is bounded by a detector 
acceptance, and the space within the acceptance is binned according to one or more bin sizes. 
Particle number n is distributed on the full momentum space. It is useful to think of transverse 
momentum p t or rapidity y t as a continuous measure distributed on angle subspace (77, 4>) 
and sampled by individual particles. Number n and transverse momentum pt or rapidity yt 
correlations can be considered both separately and in conjunction. 

An ensemble of event- wise histograms on momentum subspace x is represented schematically 
in Fig.^(first panel), with bins a, b singled out. Events can be compared with the ensemble-mean 



distribution to determine relative information, measured by fluctuations of bin contents about 
their means. The ensemble mean of event-wise pair distributions on space (x\,X2) in Fig. ^ 
(second panel) can be compared to a reference distribution consisting of cartesian products 
of single-particle mean distributions. The difference reveals correlations in the two-particle 
distribution corresponding to fluctuations in the single-particle distribution : 3_. That relation is 
the basis for the integral equation connecting fluctuations to correlations described below. 




Figure 1. Event-wise distributions on a binned primary space, corresponding two-particle 
space and possible distributions of bin contents for selected bins a and b in the primary space. 

Fig-E(third panel) sketches possible frequency distributions on number combinations (n a , rib) 
from bin pair (o, b) in the first panel. The ellipses represent half-maximum contour lines 
for gaussian-random fluctuations. The three cases correspond to correlation (solid curve), 
anticorrelation (dash-dot curve) and no correlation (dashed curve) between bins, the last 
being expected from a mixed-pair or central limit reference. The 2D frequency distribution 
is characterized by two marginal variances and a covariance. The marginal variances for bins a 
and b are given by a 2 b = (n — n) 2 b = n 2 a b — n 2 b . The covariance between those bins is given by 

a ab = ( n ~ n)a{n — n)b = n a — n a = a\ — a\^ where a\ and a\ are variances along sum and 
difference diagonals in the third panel. Pearson's correlation coefficient is a measure of relative 

covariance [I]. For bin pair (a, b) it is defined by r ab = a 2 b / J a 2 = {erf; — o"^}/! "!; + ^a) ^ 
[— 1, 1]. The numerator is the (a, b) covariance and the denominator is the geometric mean of 
the marginal a and b variances. That coefficient is our basic correlation measure. Quantities r a b, 
determined for all histogram bin pairs in the second panel, completely represent fluctuations on 
space x. Variances and covariances depend on the bin size or scale on x. The scale dependence 
of fluctuations is in turn directly related to two-particle correlations, as described in this paper. 

3. Object and reference distributions 

We distinguish between an object distribution, part of whose correlation content we wish 
to measure, and a reference distribution which contains by construction information in the 
object distribution we wish to ignore. We determine object and reference distributions for 
number n correlations and pt correlations (the latter require a slightly different treatment, as 
described below). One single-particle reference for fluctuation measurements is the ensemble- 
mean histogram on x. Variances and covariances for single bins and bin pairs measure the 
average information in the event-wise object distribution relative to the ensemble mean. 

For two-particle distributions on (x\,X2), object pair density p bj(xi,X2) is constructed from 
sibling pairs taken from same events, and reference pair density p re f{x\, 22) is constructed 
from mixed pairs taken from different but similar events. The reference could also be a 
Cartesian product of single-particle means (therefore factorizable by construction). Object 
and reference distributions are combined in several ways: 1) density of correlated pairs: 
A/o(xi,X2) = Pobj{xi,X2) — Pref(xi, X2); 2) density of correlated pairs per particle pair. 



Ap(xi, X2) / Pref{x\, X2) = Pobj(xi, X2)/ Pref(xi, X2) — 1; 3) density of correlated pairs per 
particle: Ap(x 1 ,x 2 )/^p r ef{xi,X2) = {p bj{x 1 ,x 2 ) - Pref{xi,X2)}/^Pref{xi,x 2 ). 1) and 2) 
are conventional correlation measures. 3) is unconventional, but closely related to Pearson's 
correlation coefficient described in the previous section: a relative covariance measure which 
does not depend on the absolute number of particles per se. Pearson's coefficient is ideally 
suited for testing linear superposition in the context of heavy ion collisions. 

4. Autocorrelations 

The autocorrelation concept was first introduced to time-series analysis in the form A(t) 
1 7 J f(t) f(t + t) dt, where r is the lag [5]. The concept is most useful when function 

f(t) is stationary: its correlation structure does not depend on absolute location on time. The 
information in f(t) is then fully represented by the autocorrelation distribution on r. The 
autocorrelation concept can be generalized to spatial correlations. If event-wise structure is 
randomly positioned on space x then the corresponding ensemble-average two-point distribution 
on (xx,X2) is stationary (not depending on absolute position on x). Distributions in Fig. [2 
(left two panels) of measure Ap(x\, X2)/ p re f( x i, ^2) on rj and 4> are typical of Au-Au collisions 
at RHIC 5 . Two-particle distributions can also be defined on sum and difference variables: 
p{xi,x%) — * p(x-£,xa), with xi,X2 — ► = x\ + X2,xa = x\ — X2- The data distributions in 
Fig. El exhibit stationarity — they do not depend on sum variable xs — in which case we have 
/?(££, xa) — * p( x a) t° good approximation. We then average the two-particle density over x^ 
to obtain the autocorrelation density on xa (still a 2D density, not a projection). 




Figure 2. Two-particle correlations on rj and <p f° r central Au-Au collisions at 130 GeV, 
schematic of a binned two-particle space illustrating an autocorrelation average along the k th 
diagonal and a similar averaging procedure performed directly on difference variable xa- 

The autocorrelation on difference variables can be constructed in two ways as shown in Fig. [3 
(last two panels): 1) bin space x with microbins of size e x and average the bin contents of 
space (xi,X2) along diagonals, as discussed further below, 2) bin difference variable xa in 
space (xi,X2) directly and form the corresponding pair histogram pQ. In both cases, care 
must be taken to insure that true averages are obtained and not projections by integration. 
We define autocorrelation densities p(xa), autocorrelation histograms Ak{e x ) ~ e\ p(2ke x ) and 
joint autocorrelations on two or more difference variables. If the primary distribution on x is 
truly stationary, the autocorrelation is a lossless compression of the two-particle momentum 
distribution to a lower-dimensional space. 



5. Autocorrelation density ratio for number correlations 

We now define a universal correlation measure for nuclear collisions. Returning to Pearson's 
correlation coefficient we make the following approximation 



a\ b _ {n-n) a (n-n) b (n - n) a (n - n) b _ n a n b - n a n b 



y/(n - n)l (n - n)g 

We replace the marginal number variances in the denominator of r ab by their Poisson values. The 



result is the histogram equivalent of density ratio Ap(x\, X2)/ \J p re f(%i, ^2) previously defined. 
This relative covariance can be interpreted as the number of (anti)correlated pairs per particle 
(explicit in the last expression). The third combination of object and reference distributions in 
Sec. El and the modified Pearson's normalized covariance in Eq. are equivalent. This density 
ratio is the basic correlation measure for any bin pair (a, b) on space x. We then define an 
autocorrelation in terms of density ratios for sets of bins (a, b) on (xi,X2): we combine density 
ratios as averages along diagonals k, assuming the basic distribution on (x\, X2) is approximately 
stationary. The average on index a along the k th diagonal in Fig. [2] (third panel) is 



AA k (n) _ \ {n-n) a {n-n) a+k \ _ Ap(n;ke XA ) 



Ak,ref(n) I V n a n a+k j- A ^ p ref (n; ke XA ) 



(2) 



This is the autocorrelation definition for analysis of nuclear collisions on angle space (rj, eft). This 
density ratio is an intensive correlation measure which precisely measures relative correlations 
even for excursions of object and reference densities over orders of magnitude. 

6. Extension to transverse momentum pt correlations 

We now extend the definition of the density ratio and its autocorrelation to distributions of 
transverse momentum pt on (r/,<f)). Measurement of (pt) fluctuations is described in 0. We 
treat pt as a continuous measure distributed on space x, with scalar sums of particle pt in 
histogram bins. We could write the Pearson's coefficient for transverse momentum by analogy 
with number correlations in terms of difference {pt — pt). However, the corresponding per- 
particle pt variance can be expressed as the sum of three terms: (pt — pt) 2 jn = (pt — npt) 2 /n + 
2pt (pt — npt)(n — n) jn + p\ [n — n) 2 /n, a '(^i)' variance, a pt-n covariance and a number 
variance. The three terms have forms similar to Pearson's normalized covariance, but with 
a = b defining corresponding normalized variances (also called 'scaled variances'). Each term is 
important in the overall problem of particle and pt production. We therefore want to describe 
the structure of event- wise pt distributions in terms of (pt — npt), independent of but coordinated 
with the structure of number distributions described in terms of (n — n). Pearson's correlation 
coefficient for transverse momentum correlations thus takes the form 



(n-n) a (n-n) b (p t - np t ) a {pt ~ np t ) b 

r ab 7^^= *■ 2 /- - ' W 

^/n a n b o-^ t y/n a n b 

i.e., a relative pt covariance as opposed to a relative number covariance. The geometric mean of 
marginal variances in the denominator is in this case replaced by cr? yjn a fi b , the mean of central- 
limit expectations for the pt variances. The factor crj t is however omitted in what follows to be 
consistent with the first term of the per-particle pt variance expansion above Eq. . Factors <r| t 
or p\ may be introduced as necessary in a subsequent interpretation stage. The corresponding 
Pt density-ratio autocorrelation is 



AA k i(p t -n) _ \ (p t - np t ) a b(pt - np t ) a +k,b+i I _ &p{Pt ■ n; ke VA ,le^ A ) 

= ) ,- - = ( = e ?7A e <M / W 

A k i,ref{n) { V n ab n a+k,b+i J ^ J p ref (n; ke VA , Ze^J 



That expression is formulated explicitly in terms of 2D correlations on (rj, <f>). An analogous 2D 
expression can be derived for number correlations from Eq. |2"|). 




Figure 3. Autocorrelations for 200 GeV p-p collisions determined by direct pair counting: 
number correlations (left panels) and pt correlations (right panels) and for charge-independent 
(left) and charge-dependent or net-charge (right) correlations in each case. 

Fig. |21 shows joint autocorrelations of per-particle density ratios on angle difference variables 
for number correlations (first two panels) and pt correlations (second two panels) and for 
charge-independent (like-sign + unlike-sign) pair combinations (left panel of each pair) and 
charge-dependent (like-sign — unlike-sign) pair combinations (right panel of each pair) [7j. 
Contributions from self pairs were excluded. These joint autocorrelations, determined directly 
by pair counting, provide full access to all angular correlations on (77, 4>) but are computationally 
expensive for larger event multiplicities. We can also obtain these autocorrelations by inverting 
fluctuation scale dependence, but with much less computation effort. 

7. Relating fluctuations and correlations 

Fluctuation scale dependence results from event-wise correlation structure in single-particle 
distributions which can be extracted by solving an integral equation. We first define differential 
scale-dependent fluctuation measures based on variance differences for number n and transverse 
momentum pt fluctuations, then derive the integral equation connecting fluctuations and 
correlations. Number variance difference Aa^(5x) = (n(5x) — n(Sx)) 2 /n(5x) — 1 compares 
the normalized number variance at bin scale 5x to the central-limit or Poisson expectation 1 for 
that quantity. The difference reflects correlations in the number distribution beyond that of a 
random distribution of points. The variance difference can be interpreted as the total number of 
correlated pairs (per-bin number variance a 2 ) minus the number of Poisson-correlated pairs (self 
pairs n) per particle (divided by n). The corresponding variance difference for pt fluctuations is 
Aa 2 t . n (5x) = (p t (5x) — n(5x) pt) 2 /n(5x) — u? t . Those per-particle quantities, defined in terms 
of normalized variances, are consistent in form with Pearson's normalized covariance. 

A fluctuation measurement at a single scale (STAR detector acceptance) is shown in the 
first panel of Fig. Q] j^j. The frequency histogram on random variable (pt — i ript)/(y/Wap t ) is 
compared to a central-limit reference (green dotted curve, o~p t is the ensemble-average single- 
particle variance). The pt variance difference Aa 2 t . n (5x) defined previously compares variances 
of data and reference to reveal a variance excess. Questions then arise how to interpret the 
fluctuation measurement and how to compare it to measurements made with other detectors. 

Fluctuation measurements in different bin-size or scale intervals determine different regions 
on a common distribution representing fluctuation scale dependence, as shown in the second 
panel of Fig. 9 . The variance difference from the first panel corresponds to the single point at 
the apex of the surface in the second panel. The surface is obviously structured, but what does 
the structure mean? Fluctuation scale dependence is the running integral of an autocorrelation. 
The corresponding integral equation is a linear relation between a variance difference and an 




Figure 4. (pt) fluctuations measured at the STAR detector acceptance (histogram) compared 
to a central-limit reference (dotted curve), the scale dependence of {pt} fluctuations within the 
STAR acceptance and the corresponding p t autocorrelation obtained by inversion. 

autocorrelation, including a kernel representing the binning scheme. We can express the per- 
particle variance difference on scales (5r], Sep) as 2D discrete integral 

A<J p t :n{ me Vi ne <i>) = 4 € V e <i> K mn-kl , , (5) 

k,l=i ^p re f(n;ke v ,le ( f > ) 

with kernel K mn .^i = (to — k + l/2)/m ■ (n — I + l/2)/n representing the 2D macrobin system. 
This is a Fredholm integral equation which can be inverted by standard numerical methods 
to obtain autocorrelation density ratio Apki/ \/Pref,ki as a per-particle correlation measure 
on difference variables (t^a^a) fill- The third panel of Fig. shows the autocorrelation 
corresponding to the data in the second panel, with directly interpretable structure 

8. Derivation of the integral equation 

The derivation relies on coordinating the contents of macrobins, which determine the binning 
scale for fluctuations, and microbins, which are the basis for the discrete numerical integration 
of the autocorrelation. We divide the acceptance into macrobins of varying size 5x (scale 
dependence) and microbins of fixed size e x . The scale-dependent variance is an average over all 
bins of a particular scale within the acceptance as in Fig. [5] (first panel), where the acceptance 
is Ax, the macrobin size is 5x and the macrobin number is M. The average over macrobins is 
re-expressed as an average over microbins, as shown in Fig. |S] (second panel). The average over 
microbins is re-arrange into a sum (over k) of different diagonals and an average over microbins 
on the k th diagonal. The diagonal average is the k th element of an autocorrelation histogram. 




Figure 5. Macrobins of scale 5x on a two-particle space after binning a primary space x, 
microbins of size e x relative to macrobins and two microbin schemes, one on x, the other on xa- 



The pt variance difference at scale 5x in the first line of Eq. @ is re-expressed as a 2D sum over 
microbins in the second line, with m microbins in each macrobin. Mean multiplicity n(5x) in a 
macrobin relates to the mean microbin multiplicity in that macrobin as n(5x) = mn(e x ), which 
relation is applied in the second line. The single-particle variance in the first line corresponds 
to self pairs which are excluded from subsequent pair sums. The single-particle variance term is 
consequently dropped. The third line is a rearrangement of the sum over (a, b) into a sum over 
diagonal index k and a sum over indices a, b subject to the diagonal constraint a — b = k, 
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Factors have been introduced so that the expression in square brackets is an average along the k th 
diagonal of normalized microbin covariances, with weighting factor \ln a n b jn k ~ 1. The fourth 
line identifies the square bracket as a ratio of autocorrelation histograms. Factor n k jn ~ 1 
has been absorbed into the definition of the autocorrelation ratio, but could be extracted as a 
correction factor. In the fifth line the binning system has been shifted by 1/2 bin according to 
Fig. |S] (third panel), and symmetry about the origin has been invoked (requiring an additional 
factor 2) to simplify the indexing. The histograms are finally converted to densities by including 
the microbin width. The last line, generalized to 2D on angle variables (rj,4>), is Eq. That 
integral equation provides computationally cheap 0(n) access to autocorrelations. 



9. Inversion and regularization 

Eq. ©, a discrete Fredholm integral equation, is a matrix equation of the form D = TI + N, 
with data D, image I and statistical noise N [Hlj. In principle, one could simply invert the 
matrix equation to obtain the image. However, T -1 is effectively a differentiation and acts 
therefore as a high-pass filter in the language of electrical engineering. The T" 1 N statistical 
noise term strongly dominates the image derived from a simple matrix inversion, and must be 
substantially reduced by smoothing or 'regularization' The procedure involves treating 

the image as a matrix of free values in a x 2 fit subject to Tikhonov regularization: minimize 
Xa = ||D — TI a || 2 + 1 1 L Iq, 1 1 2 , including Lagrange multiplier a which controls the role of local 
gradient operator L. The first term measures the data-integrated-image mismatch, the second 
term measures small-wavelength noise on the image. The latter is equivalent to a compensating 
low-pass filter which offsets the effect of the differentiation. The resulting image is represented 
by Iq, = T^ 1 (D — N„) which estimates true image I. 

Regularization is illustrated in a power-spectrum context by Fig. El (first panel): a tradeoff 
between information loss and noise suppression. The basis for choosing the optimum a is 
illustrated in the second panel. ||D — TIq,|| 2 (dots) is signal loss and ||LI a || 2 (triangles) is 
residual noise. As noted, a controls a compensating low-pass filter: small values retain all signal 
and a large amount of noise from the differentiation, larger values reduce noise, and finally 
distort the signal by over smoothing. The optimum value is determined as in Fig. (second 



panel). This example illustrates that there are clear criteria for choosing an optimum a so that 
negligible information is lost from the image while statistical noise is greatly attenuated. 
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Figure 6. Sketch of Fourier power spectrum, distribution of two terms in x 2 on Lagrange 
multiplier a, autocorrelation from fluctuation inversion and corresponding smoothing error. 



Statistical and systematic errors are determined by looping through integration and 
differentiation twice in the sequence D — ► l a — > D Q — » Iq, q , including inversion, forward 
integration and second inversion. Difference D — D a estimates statistical error on the data 
and may itself be inverted to determine residual statistical error on image I Q . Difference 
I» — lao estimates the smoothing error. Fig. (right two panels) shows image I Q and 
corresponding smoothing 6iror I a — for data from the Hijing JVIonte Carlo |12j . dominated, 
by jet correlations [S]. 




Figure 7. Autocorrelations from the Pythia Monte Carlo for charge-independent (left panels) 
and charge-dependent (right panels) correlations and direct (left) vs inversion (right) methods. 



A precision comparison of autocorrelations obtained directly by pair counting and indirectly 
by fluctuation inversion is presented in Fig. Q The data were obtained from the Pythia 
Monte Carlo 13 j. The four panels in sequence are charge-independent direct and inverted 
autocorrelations and charge-dependent direct and inverted autocorrelations. The agreement 
between methods is excellent. 



10. Correlation types 

Fig. El illustrates types of correlation measurements. The left panels show an autocorrelation 
obtained by {pt} fluctuation inversion (first panel) and a sketch of corresponding two-bin 
correlations (second panel) for an arbitrary bin separation. In this type of correlation, structure 
may occur in most or all events, but with random position on (77, ^). The autocorrelation of a 
positive-definite measure (such as n or pt) must be positive-definite at the origin (a variance). 
Elsewhere (covariances), positive autocorrelation bins indicate correlation (solid curve in second 
panel), negative bins indicate aniicorrelation (dash-dot curve), as in the magenta areas adjacent 
to the positive same-side peak in the first panel. The dashed circle represents the reference. 

The right panels illustrate a case where localized structure with approximately fixed position 
is present in a minority of events. The third panel shows correlations on transverse rapidity 



Figure 8. pt autocorrelation with regions of positive and negative covariance, distributions 
of bin-pair contents illustrating corresponding correlation (solid) and anticorrelation (dash-dot) 
trends, distribution on (yti,Vt2) illustrating soft and hard components, and distributions of bin- 
pair elements illustrating the role of rare hard events in producing positive-definite covariance. 

yt (that distribution is not an autocorrelation). In the fourth panel distributions are sketched 
for two event classes. The dashed curve represents common (e.g., soft) events, the solid circle 
represents exceptional events (e.g., hard events which contain a detectable parton scatter) from 
an ensemble of p-p collisions. Positive covariances result from the exceptional events which 
have an additional multiplicity contribution localized on the yt space and (in contrast to the 
autocorrelation example at left) occuring at a nearly fixed position over the event ensemble. 
The normalized covariance density comparing sibling and mixed pairs in the third panel reveals 
the contribution from exceptional events (parton fragments). 

11. Autocorrelations and conditional distributions 

Conventional studies of jet correlations in A-A collisions employ a leading-particle analysis 
(invoked when full jet reconstruction is not possible) [Tl]. The goal is to estimate a parton 
momentum by that of the highest-pf (above some threshold) particle in a collision — the leading or 
trigger particle. The analysis utilizes two or three conditional distributions as illustrated in Fig. El 
(first two panels) , where for the sake of comparison transverse momentum pt has been replaced by 
transverse rapidity yt. The condition on (yti,yt2) is a rectangle representing asymmetric trigger- 
and associated-particle yt conditions. Trigger region (f2t tl , ^t t2 ) is by construction displaced from 
the sum diagonal. Angular correlations are also defined as conditional distributions relative to 
the trigger-particle position in single-particle angle space (77, (/>). The angular correlation is 
plotted (using pseudorapidity as an example) on conditional angle difference A 77 = tj — fjtriggeri 
which represents an event-wise shift of the single-particle angle origin. 




Figure 9. Conditions defined on (yn,yt2) and (771,772) for a leading-particle (trigger) analysis 
compared to a condition on (yti:,ytA) and autocorrelation on va invoking no leading particle. 

Alternatively, one can abandon attempts to estimate parton pt per se and use a technique 
involving no conditions on momentum, or a symmetric condition on (yti,yt2) with no conditions 



on angle variables, as illustrated in Fig. (second two panels). The fragment distribution on 
{yn-)Ut2) can be studied in its own right. Cut conditions on (yti,yt2) can be defined to study 
corresponding changes in jet angular morphology on (rj, <f>). Angular autocorrelations invoking 
no trigger condition and defined on symmetric difference variables such as i]A = Vi ~ V2 access 
a minimum-bias parton population. Difference variable t/a spans the diagonal axis of the 2D 
(Vii 7 ]?) space. Variables Arj and 77A are numerically different, support different distributions and 
should not be confused (the same comment applies to Acj) and 4>a)- The correlation measures 
used in the two cases may be compared: 

Ap(n;r? A ,<M) f d 2 N pair = p(n; Arj, A<j)) 

■J p re f(n; r/A 5 <^a) ^ N **B dAr] dA(/) ~ ^tn g d P* P( n '> P^d) ' 

The latter depends directly on leading- or trigger-particle acceptance tltrigi an d a reference must 
be provided a posteriori by defining a model function to describe the background. 

12. Summary 

We have derived an integral equation which connects fluctuation scale dependence to 
corresponding autocorrelations. Inversion of the integral equation reveals autocorrelations which 
are equivalent to those from pair counting. Fluctuations are thereby interpretable in terms 
of underling two-particle correlations. Autocorrelations are complementary to leading-particle 
techniques for analysis of jet correlations. Definition of normalized variances, covariances and 
variance differences is tightly constrained by a number of considerations, including linearity, 
relating fluctuations to correlations, minimizing statistical bias and insuring that correlation 
structure is maximally accessible and interpretable over a broad range of contexts. 
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